Topological edge states in the one-dimensional super-lattice Bose-Hubbard model 
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We analyze the ground state of interacting ultra-cold bosonic atoms in a one-dimensional (ID) super-lattice 
potential with alternating tunneling rates ti and t2- A topological order parameter is introduced which is quan- 
tized in units of tt for the Mott insulating (MI) phases. A step in the effective confining potential created e.g. by 
a second heavy atom species can lead to an interface between two MI regions with filling n = 1 and n = 1/2. 
Depending on the ratio ti/t2 the n = 1/2 MI phase is topologically non-trivial which results in localized, 
protected many-body edge states. Density-matrix renormalization group (DMRG) simulations show that the 
edge states manifest themselves either in localized density minima or localized density maxima at the interface, 
which can easily be detected. Shape and energy of the edge states as well as conditions for their occupation are 
determined analytically in the strong coupling limit and in general by DMRG simulations. 

PACS numbers: 05.30.Fk, 03.75.Hh, 73.21. Cd 
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Topological phases have become an intensively studied 
subject in many fields of physics. Key features of condensed- 
matter systems such as topological insulators Jll-Ill or su- 
perconductors 10] as well as quantum Hall systems IH-lsl 
have been related to the existence of robust edge states at 
interfaces between phases with different topological charac- 
ter. Currently there is a growing activity in condensed-matter 
physics aiming at creating and understanding novel topolog- 
ical phases. On the other hand ultra-cold atomic gases have 
developed into an ideal experimental testingground for con- 
cepts of solid-state and many-body physics ^ and they could 
become important for studying topological effects lUOl - llSll . 
One of the simplest models possessing non-trivial topologi- 
cal properties is the Su-Schrieffer-Heeger (SSH) model Ill4ll . 
which can be realized by ultra-cold fermions in a ID tight- 
binding super-lattice (SL) potential with alternating hopping 
amplitudes. Its topological properties on a single-particle 
level have recently been explored both theoretically [fisl [l6ll 
and experimentally lUsll . Here we show that in the case of 
bosons, MI phases with filling n = 1/2 can be non-trivial 
topological insulators. In this case, introducing a localized po- 
tential step, realized e.g. by the admixture of a second, heavy 
atomic species, allows to create an interface between gapped 
MI phases with different topological invariants. We show that 
at the interface, a localized and robust many-body edge state 
emerges which manifests itself in density minima or maxima 
that can easily be observed with techniques developed in re- 
cent years lUTMlOll . The edge state, which is an indicator of 
the topology of the bulk phases, is shown to survive until the 
insulating phases melt due to tunneling. 

The starting point of the discussion is the ID SL-Bose- 
Hubbard model (SL-BHM) in the grand-canonical ensemble 
described by the Hamiltonian K ^ H ~ 
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FIG. 1 : (Color online) (a) Phase-diagram for the SL-BHM with t2 = 
0.2ti obtained by DMRG. One recognizes the presence of MI phases 
with integer and half-integer filling, (b) Different dimerizations I and 
II of SL potential corresponding to Zak phases i^' = and 



the bosons at lattice site j, and fij 



Particles can tunnel 



where ctj and aJ- are the annihilation and creation operators of 



from a lattice site to a neighboring one with alternating hop- 
ping amplitudes ti and t2 and there is an on-site interaction U. 
Ej describes a potential which can also include a slowly vary- 
ing harmonic confinement, and ^ is the chemical potential. 
Eq. ^ corresponds to the SSH model in the case of fermions 
(where the interaction term is absent) or in the equivalent case 
of hard-core bosons {U — > oo). 

A generic ground-state phase diagram of the SL-BHM is 
shown in FigH] which is taken from Ref . ll20ll . Besides MI 
phases with integer filling, known from the simple BHM, it 
shows loophole-type insulating regions with half integer fill- 
ing for ti > t2 11211 122I1 . These phases shrink in size when 
ti decreases, and vanish at ti = t2, i.e. in the case of simple 
BHM. They reappear when ti < t2 and the point ti = t2 
marks a topological phase transition. 

In the case of non-interacting fermions, the topology of the 
band structure is determined by its Zak phase ll23ll (or wind- 
ing number), u ~ i J^^ dk {u{k) \ dk \u{k) ), where \u{k)) are 



2 



y-1 IT! 


n 74- 

U. 1 4: 


n 9^ 






\vi /tt 

2 

\ui /tt — 1 


3.1 • 10"^=* 
7.9 ■ 10"^^ 


6.1 • 10"^^ 
3.0 • 10"^^ 


1.0 ■ 10"^^ 
8.9 ■ 10~" 


4.4 • 10"^^ 
1.4 ■ lO"^'^ 


U 5 lit 


-0.27 


0.54 


-0.12 


0.85 


V h I'K 
12 


-0.042 


0.13 


-0.46 


0.0082 



TABLE I: Numerical values of the many-body winding number v, 
for a n = 1/2 MI of A'' = 5 particles {y\, upper block) and for 
N — 5 particles on 12 sites (i^_5_, lower block) corresponding to a 

compressible phase. In the MI case, for both 1 2 /i i = 2 (upper lines) 
and t2/ti — 1/2 (lower lines) v is quantized in units of tt within the 
numerical precision. 



the single-particle Bloch functions and k is the lattice quasi- 
momentum. While for a general ID band structure this phase 
may take arbitrary values, it is integer quantized (in units of tt) 
112411 when the sub-lattice symmetry is present. In the case of 
the SSH model the winding numbers of the upper and lower 
band are equal but opposite 1^ = (±7r) ifssll for dimerization 
I (II), see Fig. □ 

For interacting systems there is no conserved lattice quasi- 
momentum k and we must construct a many-body generaliza- 
tion of the winding number. Like the Chem number llsoll it can 
be defined via generalized boundary conditions, ip{xj + L) = 
e^^/'"'ip{xj) for all particle coordinates j ~ 1, N and sys- 
tem size L. They correspond to a magnetic flux 6 threading 
the system. When this flux is adiabatically varied, the many- 
body wavefunction \ '^{9) picks up a Berry phase [i3 lil 



de {ij{e)\de\i^{e)). 



(2) 



This topological order parameter can easily be calculate in the 
hard-core limit U — ?> 00. In this case the n = 1/2 MI phase is 

weU described by 1*1/2) = j^^ dx wii^x ~ Rj)^"^ {x) |0), 
where wi denotes the Wannier function of the lower SSH band 
and the lattice vectors Rj are Wannier centers. A similar ex- 
pression holds for the n = 1 MI phase, where also the higher 
SSH band is filled. We find that the Mis with integer fin- 
ing are topologically trivial with v = Q and those with half- 
integer filling can take the values = Q for dimerization I 
and for dimerization II. 

Most importantly numerical exact diagonalization of small 
systems shows that the topological invariant stays strictly 
quantized even for finite U as long as the particle-hole gap 
remains finite (see table |I]i. Thus we expect that the non- trivial 
topology of the SSH bands for free fermions carries over to 
bosons with finite interactions. This is our motivation to study 
edge states between topologically distinct MI phases in the 
SL-BHM as indicators for a quantized topological invariant. 

An interface can be created in ultra-cold lattice gases as fol- 
lows: An increase of the potential energy in the Hamilto- 
nian ([T]i corresponds to a shift in the chemical potential. Thus 
an increase of €j for a number of consecutive lattice sites by 
Ae can lead to an interface between two MI phases with inte- 
ger (e.g. n = 1) and half-integer filling (e.g. n — 1/2). This 
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FIG. 2: (Color online) Ground-state density distribution of the SL- 
BHM with harmonic trap and additional potential step between sites 
j = — 6 and j = 5 leading to interfaces between n = 1/2 (in 



the center) and n = 1 MI regions, e 



{j + 0.5y, with 



ui/U = 0.001. Results are obtained by DMRG simulations for 
fi/U = 0.55. (a) topological trivial 71 = 1/2 MI phase with 
ti/U = 0.04, t2/U = 0.2 and Ae/U = 0.6. (b) topological 
non trivial n = 1/2 phase with ti/U = 0.2, t2/U = 0.04 and 
Ae/U = 0.6(blue squares), 0.7(red stars). Right panel illustrates 
interface in topologically trivial case (c) and in nontrivial case with 
occupied (d) and unoccupied interface (e). 



can be realized if n and Ae are chosen such that 



fil < fi < fi\ 



1/2 



Ae < ^ < /i 



1/2 



Ae 
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Here f/J^ and denote the upper (+) and lower (— ) bound- 
ary of the insulating regions in the phase diagram of Fig{T] As 
shown e.g. in Ref.fSJl for the case of bose-fermi mixtures, 
an effective potential step can be created by an admixture of a 
second atomic species, e.g. fermions, with very small hopping 
rates. Under appropriate conditions (see ||32|) a repulsive, lo- 
cal interaction V between the two species results in an effec- 
tive attraction of the fermions. With a weak harmonic con- 
finement the fermions form a connected cluster at the center 
of the trap with unity filling and sharp boundaries if the result- 
ing binding energy exceeds their hopping. This results in an 
increase of the potential energy of the bosons Ae = V which 
extends over all sites of the fermion cluster If Ae and fi fulfill 
(|3]l and (|4]i, two interfaces between an = 1/2 MI at the cen- 
ter of the trap and surrounding n = 1 MI regions are created. 
Depending on the location of the interfaces relative to the sub- 
lattices the winding number either stays, Ai/ = 



'1/2 



or jumps. 
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Fig. 12] shows the density distribution in a weak har- 
monic trap with an additional potential step Ae calculated 
by density-matrix renormalization group algorithm (DMRG) 
II25I - I27II . One clearly recognizes interfaces between a central 
ri = 1/2 MI and surrounding n = 1 MI regions. Since the 
number of heavy particles was taken to be even, both inter- 



faces are characterized by the same change | Aj/| in the wind- 
ing number. The upper plot shows the case Av = 0, the lower 
one /S.V — TT. In the first case there is a simple step in the den- 
sity and no additional structure at the edge. The same holds at 
an interface between any two MI phases with integer fillings 
irrespective of the dimerization. In the second case, however, 
one sees pronounced dips or peaks in the average density. We 
will argue now that these are clear signatures of the edge state 
which must be present at any interface between topologically 
distinct bulk phases. 

To understand the origin of the topological edge state, we 
note that in the hard-core limit [/ — > oo our situation corre- 
sponds to free fermions and an interface between a half filled 
SSH band (n = 1/2 side) and completely filled SSH band 
(n = 1 side). With an open boundary and for dimerization II 
(see FiglT})) a state with energy E ~ Q appears between the 
two half bands which are separated by a finite gap for ti 7^ t2, 
while for dimerization I there is no such state. If only the 
lower SSH half-band is filled the mid-gap state leads either 
to a peak in the density distribution at the boundary if it is 
occupied by a fermion or to a dip if it is not occupied. On 
the other hand for complete filling no such structure exists, 
as is the case for dimerization I. This is the essence of the 
edge state counting formula lfl6ll which states that at an open 
boundary of a topologically nontrivial phase with v = t: ex- 
actly one localized edge state exists, which can be occupied 
by a fermion or not. If the open boundary is replaced by an 
interface between a half filled and a fully filled band, the finite 
tunneling ti leads to a renormalization of the energies of the 
two mid-gap states, which remain however within the corre- 
sponding gaps. In the following we discuss what happens if 
the interaction U becomes finite which leads to a mixing of 
single-particle states. 

We start by discussing the properties of the edge states in 
the limit of large interactions U ^ ti and small hopping 
t2 ^ ti. This allows a perturbation expansion in t2 and 
1/y (" cell strong coupling perturbative expansion", CSCPE 
1128112911 ). We note that the same analysis applies for an inter- 
face between any 2i±i and ^ Mis, with m being an integer 
In this case the low lying bands can be assumed inert and only 
lead to a renormalized hopping due to Bose enhancement. We 
discuss the grand-canonical ground state and consider states 
-01 and ■00 with and without an additional particle at the in- 
terface, see Figs|2}l and e. The critical chemical potentials 
defining the borders of the MI phases in the bulk in FiglUcan 
easily be calculated in CSCPE. We find up to terms of order 
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FIG. 3: (Color online) The red line, = E^i^ — Ei'^\ separates 
parameter regions with an occupied edge state (upper part) or unoc- 
cupied edge state (lower part). When it is within the shaded region, it 
also indicates stability of the edge states towards hybridization with 
the bulk. Parameters are ti/U — 0.2 and *i/*2 = 5. 



To understand whether the edge state is actually bound to the 
interface or hybridizes with the bulk, we perform a stability 
analysis. We therefore investigate under what conditions the 
many-body state lies within the gap to collective excita- 
tions on the ground state ipQ. In this case both many -body 
states i/'o,i constitute a localized particle or hole at the edge 
and we call the edge state stable. By tuning yu both states can 
be made the grand-canonical ground state. Hybridization can 
occur if a particle or a hole excitation is created in the bulk of 
either of the two Mis. Reducing (increasing) the number of 
particles at the MI interface by one and simultaneously creat- 
ing a particle (hole) in the bulk costs a finite energy if 



Here denote the energies of i/jq^i and CSCPE yields 



' Ae < Ei'^ - eP < u]!' + Ae 



(5) 
(6) 
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+ 4*2 + C»(*^/C/,l/[/3) (7) 



From Figl3]one can read off the potential heigths Ae for 
which non of the states ?/)o,i hybridizes with the bulk: The red 
line has to be within the shaded area. Moreover, ii'i) is the 
grand-canonical ground state for a chemical potential fjL < 
(fi > ^c) with = eP 



eP . For yu > /ic, i.e. above the 



1/2 + 4- 

lij =*2-*i. 



red line in Figl3] the edge state is occupied, while below it is 
not. The first case (ipo) can be detected by measurement of a 
density less than 1/2 on the n = 1/2 MI side of the interface, 
the second (ipi) by measurement of a density larger than 1 on 
the n = 1 MI side. 

To verify these results we performed numerical simulations 
for a step potential ej/U = AeO(j — jstep + 0.5), where O 
is the Heaviside step function. In Fig|4^ we show DMRG 
results for the density as a function of /i. For iJ,/U = 0.45 
the system is inside the stability region of Figl3]above the red 
line. One clearly recognizes a well localized dip in the den- 
sity. /i/C/ — 0.50 corresponds to an occupied edge inside 
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FIG. 4: (Color online) (a) Density distribution at potential step 
Ae/[/ = 0.6at jstcp = 6forti/[/ = 0.2, ti/t2 = 5 and increasing 
chemical potential /i/(7. (b) Density {n\ °'^) at the edge of the n = 1 
MI region as function of fi/U showing the occupation at the edge if 
fi > He = 0.47 (dashed red line). 

the stability region. Here a clearly pronounced density peak 
appears. When /i is chosen such that the system is outside 
the region of the n = 1/2 Mott insulator (fi/U = 0.37 and 
fi/U = 0.65) the density dip on the n = 1/2 side starts to 
vanish while the peak on the ?? = 1 MI side remains. Fig|4j5 
shows the density (nf^^) at the edge of the n = 1 MI side as 
function of i-i/U. As soon as jj, exceeds (red dashed line), 
there is a clear jump indicating the transition from occupied 
to unoccupied edge. 

A particular feature of the edge state in the SL-BHM, not 
present in the hard-core limit, is the peak of the local density 
on the border of the ?? ~ 1 Ml region above unity. Within 
CSCPE we can calculate this density, which yields in zeroth 
order of t2 

Fig|5t shows a comparison of (nf'^'') obtained by DMRG 
with eq.(|8ll (green, dashed line). Also shown is the CSCPE 
result in second order in t2- The corresponding analytic ex- 
pression can easily be obtained but is too involved to be given 
here. One recognizes very good agreement even for rather 
large tunneling rates just before the n = 1 Ml starts to melt. 
For ti/U —5- 0, (n^'^^'^) approaches the hard-core value, while 
finite tunneling leads into {fil'^^'^) > 1. Beyond some opti- 
mum value of ti/U the edge state starts to delocalize as the 
particle-hole gap of the insulator becomes smaller and {h^'^^^) 
decreases. Finally Fig|5j5 shows density distributions for the 
occupied and unoccupied edge for ti/U = 0.5, i.e. close to 
the melting point of the n = I Ml phase. One clearly sees that 
the edge state is still stable and localized. 

In summary we have discussed topological properties of the 
ID superlattice Bose-Hubbard model with alternating hop- 
ping rates ti and ^2- In the limit of infinite interaction U 
this model corresponds to the SSH model for free fermions, 
which is known to possess topologically non-trivial insulating 
phases for ti ^ t^- We introduced a many -body generaliza- 
tion of the Zak phase as topological order parameter Based 
on small system exact diagonalization we argued that also for 
the bosonic model with finite interactions this order parame- 
ter is strictly quantized for the Mott insulating phases. To ver- 
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FIG. 5: (Color online) (a) Density (^1'°^) at the edge of the n = I 
MI region as function of ti/f/ extracted from DMRG (blue squares) 
compared to analytic results, eq.ljS) (green dashed line), and includ- 
ing second order terms in t2 (red full line). Other parameters are 
fi/U = 0.96 - l.Qti/U, Ae/U = 0.88 - OMi/U and h/tz = 5. 
(b) Density distribution at interface between topologically trivial 
n — 1 MI phase and topologically non-trivial n = 1/2 MI phase 
obtained by DMRG for 1 1 /f/ = 0.5, i 1 /t2 = 5 and js,ep = 16. The 
inset shows the situation with ti and t2 exchanged, corresponding to 
an interface between two topologically equivalent phases. 



ify their topological nature we analyzed edge states between 
a topologically nontrivial MI with filling n = 1/2 and a MI 
with filling n = 1 using DMRG and analytic perturbative cal- 
culations. At the interface a stable and well localized edge 
state emerges, whose properties were studied in detail both 
analytically and numerically. Our numerical simulations ver- 
ified that the many-body system supports exactly one local- 
ized, excited state within the bulk gap of the Mis. We argue 
that this is a result of the non-trivial topology of the n ~ 1/2 
MI phase and notice that the edge state counting formula 1II6I1 
carries over to the many-body winding number In a grand- 
canonical setting the edge can have two different occupation 
numbers or 1 in the interesting range of chemical potentials. 
The first case is characterized by a density dip at the edge of 
the 71 = 1/2 MI region below 1/2. For finite values of V the 
second case results in a density peak at the edge of the n = 1 
MI region with a local density exceeding 1 with a maximum 
value for intermediate values of tijlJ . These features allow a 
simple detection of the edge states and thus a verification of 
the different topological nature of the MI phases in cold-atom 
experiments. The required potential step can be realized by an 
admixture of a second heavy atom species in the presence of 
a weak harmonic confinement potential. 
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